Hahn echo and criticality in spin-chain systems 
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We establish a relation between Hahn spin-echo of a spin-| particle and quantum phase transitions 
in a spin-chain, which couples to the particle. The Hahn echo is calculated and discussed at zero as 
well as at finite temperatures. On the example of XY model, we show that the critical points of the 
chain are marked by the extremal values in the Hahn echo, and can influence the Hahn echo in finite 
temperatures. An explanation for the relation between the echo and criticality is also presented. 
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Quantum phase transitions l](QPTs) have attracted 
enormous attention within various fields of physics in 
the past decade. They exist on all length scales, from 
microscopic to macroscopic. Because QPTs, which de- 
scribe transitions between quantitatively distinct phases, 
are driven solely by quantum fluctuations, they provide 
valuable information about the ground state and nearby 
excited states of quantum many-body systems. The ob- 
servation of quantum criticality depends eventually on 
the experimentally available temperature, then it is nat- 
ural to ask how high in temperature can the effects of 
quantum criticality persist? Do quantum critical points 
shed light on quantum mechanics of macroscopic sys- 
tems, for instance providing a deeper understanding of 
decoherence? In this paper we answer these questions by 
examining exact solutions of the XY spin-chain model, 
calculating the Hahn echo of a spin- 1/2 particle coupled 
to the critical spin-chain. 

The Hahn echo was first introduced by HahnQ to ob- 
serve and measure directly transverse relaxation time 
T%, i.e., the dephasing time. It differentiates from the 
Loschmidt echo in that the latter measures the sen- 
sitivity of quantum system dynamics to perturbations 
in the Hamiltonian. For a certain regime of parame- 
ters, the Loschmidt echo decays exponentially with a 
rate given by the Lyapunov exponent of the underly- 
ing classically chaotic system. Recently, a huge inter- 
est was attracted in the attempt of characterizing QPTs 
in terms of entanglement, by analyzing extremal points, 
scaling and awmptotic behavior in various entanglement 
measures 

[1 II U H The relation between Berry's 
phases and quantum critical poin ts was also established 
recently in the XY mo del [sllii. I lOf. In this paper, we shall 
show how critical points can be reflected in the Hahn 
spin-echo, and what is the finite temperature effect on 
the Hahn spin-echo. 

Consider a spin-1/2 particle coupled to a spin-chain 
described by the one-dimensional XY model, the Hamil- 
tonian of such a system may be given by 



H — H s + H c + Hi 



(1) 



where 



JV 



H c = -2^((l + 7 ) s f s f +1 + (l- 7 ) s f s f +1 +A S f), 
1=1 

N 

4<? (2) 



H 
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Here s denotes spin operator of the system particle which 
couples to the chain spins si (1 = 1,..., N) located at the 
lattice site I. The spins in the chain are coupled to the 
system particle through a constant g. The Hahn echo 
experiments consists in preparing the system spin in the 
initial state \y s ) = (| T) + A J.))/v2> and then allowing 
free evolution for time r. A 7r-pulse described by the 
Pauli operator a x is then applied to the system spin, and 
after free evolution for one more interval r an echo is 
observed, which provides a direct measurement of single 
spin coherence. We would like to notice that the free 
evolution here means no additional driving fields exist, 
the coupling between the system and the spin-chain is 
always there. 

We now follow the calculation|ll| to derive an exact 
expression for the Hahn echo decay due to the system- 
chain couplings in Eq.Q. The density matrix for the 
whole system (system and the spin-chain) which will be 
used to calculate the Hahn echo is given by 



p{r) = U{t) Pq U\t), 
where U (r) denotes the evolution operatorfTl] 



U(t) = u(T)a x u[r),u(T) = e 



-iHr 



(3) 



(4) 



and po is taken to be the initial state of the whole system 

po = \y s )(y s \®Pc(0) (5) 

with p c (0) denoting the initial state for the spin-chain. 
The Hahn spin echo envelope is then given by 



v e (t) = 2\Tr{(s x + t S y)p(r)}\. 



(6) 



In order to get an explicit expression for the Hahn echo 
envelope Eq.lJBJ, we first write u(t) in basis {| |), | I)} 
(the eigenstates of <r z ), by noting that [H s ,Hi\ = 0. This 
leads to 



ps' 



i=T,i 



2 



with Uj(r) satisfying, 

d 

N 

H, = -2]T((1 + 7 ) s f S f +1 + (1 - 7 ) S f S f +1 + Xjsf), (8) 

2 = 1 

where Xj = X±g, + and — correspond to J. and T, respec- 
tively. The free energy /x of the system which contributes 
only energy shifts to Hj would not affect the Hahn echo 
and has been omitted hereafter. For the system initially 
in state \j) (j =t, J.), the dynamics and statistical prop- 
erties of the spin-chain would be govern by H j , it takes 
the same form as H c but with perturbed field strengths 
Xj. This perturbation to the spin-chain regardless of 
how small it is can be reflected in the Loschmidt spin 
echo decay 0], m particular at critical points. What 
behind the decay is the orthogonalization between two 
ground states obtained for two different values of external 
parameters [13|. The Hamiltonian Hj can be diagonalized 
by a standard procedure to be 

k 

which can be summarized in the following three steps. (1) 
The Wigner- Jordan transformation, which converts the 
spin operators into fermionic operators via the relation 
a i = (Ilm<2 a m)( a i + *°f where 07 is the Pauli ma- 
trix of the spin at site l; (2)The Fourier transformation, 
dk = a/exp(— i27r/fc/iV); And (3) the Bogoliubov 

transformation, which defines the fermionic operators, 
Vj,k = d k cos id'_ k sin-^-, (10) 

where the mixing angle Qj.k was defined by cos 9j ik = 

£j,k/uj,k, with oj jtk = \Je 2 j k + j 2 sin 2 and s jtk = 

(cos^-Xj), k = -iV/2,-JV/2 + l,...,iV/2-l.Todiag- 
onalize the spin chain Hamiltonian, the periodic bound- 
ary condition was used in this paper. The boundary 
term H b = -(a\a N + a^ai + j(a N a\ + h.c.))(P + 1) 

with P = exp(iTT^2j =1 ajjaj) would vanish when N/2 is 
odd. Since the paper aims at finding the link between 



the Hahn echo and the critical points, we will chose N/2 
odd to simplify the boundary effects. This treatment is 
available in the limit N — > 00, where the boundary effects 
are negligible. It is easy to show that [jy^fc, rjj^] ^ when 
j i, i.e., the modes rji t k and rjj^ k do not commute(this is 
not the case for some special parameters discussed later 
on). This would result in the Hahn echo decay as you 
will see. With these results, the evolution operator U(t) 
can be reduced to 

U(t) = u t (tH(t)| T)U I + «i(r)« T (r)| |><| |, (11) 

with u.j{t) = e^E^.^L^-^ = n fc «i,fc(r), j =T 
,|, and Uj,fc(r) = e'^^t"^"'^ . After a simple al- 
gebra, we arrive at 

ve(t) = |Tr c (u{u|u;u T p c (0))|, (12) 

where the trace is taken over the spin-chain. Ea. (|12f> can 
be simplified by noting that (rij.k — v] k 1 !]^) 

[n h k,n Lk ] = -sin(0 Tlfe - 9 Lk )(v-kVk - vW- k )> ( 13 ) 

and consequently, 

u\ k ui,k = u Lk u\ k + X k , (14) 

where X k = (1 - e M T, fe*)(l _ e~ lu} ^ kt )[n hk , n ltk \. Here 
i] k = d k co s % - i(iL fc sin^, and 8 k = 6j ik \\ j= \. Substi- 
tuting Ea. 1)1 into ve(t), we get 

ve(t) = I] I 1 + ^r c [u\ k X k u hk pM\\- (15) 

k 

The explicit expression for Ea. p5|l can be obtained by 
choosing a specific initial state of the chain. We shall 
consider two initial states in this paper, (1) p c (0) is taken 
to be the ground state of H c , (2) p c (0) is chosen to be 
a thermal state for the spin-chain. The ground state of 
H c follows by the same steps summarized above. It is 
defined as the state to be annihilated by each operator 
rjk, namely ?7fc|<?(7, A)) = 0. After a few manipulations we 
obtain the Hahn echo envelope at zero temperature, 



v E (t) = I 1 + 4 sin ( T,fc - h,k) sm(0 fe - T ,*)|1 - e^' fc T (1 - e^*' - sin 2 k - Lk \l - e ^f ) 

k 

- 1 sin^ - 6 Uk ) sin(0 fe - e Lk )\\ - e*^*<| 3 (l - e*T.»* - sin 2 9 ±Z±A \i _ e *"t.*f)|. (i 6 ) 

I 

With the above expressions, we now turn to study the is exactly solvable and still present a rich structure, it 
Hahn echo at zero temperature. Since the XY model offers a benchmark to test the properties of Hahn echo 
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FIG. 1: (color online) Hahn echo of the spin-^ particle vs. 
time r and the anisotropy parameter 7. The spin-i particle 
was coupled to a spin-chain described by the XY model. The 
parameters chosen are N = 246 sites, g = 0.3 and (a)A = 2, 
(b) A = 1.5, (c)A = 1, and (d)A = 0.5. 

in the proximity of a quantum phase transition. For the 
XY model one can identify the critical points by finding 
the regions where the energy gap ujk vanishes. Indeed, 
there are two regions in the A, 7 space that are critical. 
Namely, 7 = for — 1 < A < 1, and A = ±1 for all 7. 
We first focus on the criticality in the XX model. The 
XX model that corresponds to 7 = has a criticality 
regime along the lines between A = 1 and A = — 1. The 
critical points can be read out from the Hahn echo as 
shown in figures and [21 Figure 1 shows the Hahn echo 
as a function of time r and the anisotropy parameter 7. 
Clearly, the Hahn echo takes a sharp change in the limit 
7^0, this results can be understood by considering the 
value of 9j.k and 9k, which take or ir depending on the 
sign of cos(27rfc/7V) — Xj and cos(2irk/N) — X, respectively. 
In either case, sin(9k — Oj^) = sin(0|^ — O^k) = 0, this 
leads to ve(t) — 1. Physically, when 7 = 0, the particle 
number operators Tif k and commute, which implies 
that the perturbation from the system to the spin-chain 
does not excite the spin-chain, then the Hahn echo which 
characterizes the dephasing of the system remains unit. 
Figure 2 shows the Hahn echo ve(t) in the vicinity of 
critical points 7 — > and A = ±1. A sharp change among 
the line of A = ±1 appears clearly. 

We would like to notice that the Hahn echo ve(t) at 
critical points of 7 = and A = ±1 does not depend on 
the chain-system coupling constant g, but in the vicinity 
of 7 = 0, it does. This was shown in figure 3, where 
we plotted the Hahn echo as a function of A and g with 
7 = 0.001 (close to zero). As expected, the critical points 
have been shifted linearly by the coupling constant g. 
The white area in figure 3 corresponds to ve{t) = 1. 
In the region of g > 2 and —1 < A < 1, ve(t) always 
equal to 1 . This can be understood by examining the the 
definition of 9j k and 9k- In this region, Xj = X ± g > 1, 
leading to 9j t k = 9k for any k in the limit 7 — > 0. This 
results in Ve{t) = 1, which is a direct followup of Ea. (|16ll . 

Now we turn to study the criticality in the transverse 
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FIG. 2: (color online)Hahn echo as a function of time r and 
A. The figure was plotted for N = 246 sites, g = 0.1 and 
7 = 0.001. 




FIG. 3: (color online)This figure was plotted to show the 
dependence of the critical points on the system-chain coupling 
constant g. Time r = 50, N = 246 sites, and 7 = 0.001 were 
chosen for this figure. 

Ising model(7 = 1 in the XY model). The ground state 
structure of this model change dramatically as the pa- 
rameter 7 is varied. We first summarize the ground states 
of this model in the limits of |A| — > 00, |A| = 1 and A = 0. 
The ground state of the spin-chain approaches a prod- 
uct of spins pointing the positive/negative z direction 
in the |A| — > 00 limit, whereas the ground state in the 
limit A = is doubly degenerate under the global spin 
flip by Yli—x erf- At |A| = 1, a fundamental transition in 
the ground state occurs, the symmetry under the global 
spin flip breaks at this point and the chain develops a 
nonzero magnetization (a x ) ^ which increases with A 
growing. The above mentioned properties of the ground 
state are reflected in the Hahn echo as shown in figure 
4. In the limit |A| — > 00, 9j^ = 9k = 7r/(— 7r), this re- 
sults in ve(t) = 1. In fact as figure 4-(a) shows, when 
\M > 4, ve{t) approaches 1 very well. With |A| — > 1, the 
Hahn echo ve{t) tends to zero, this can be interpreted as 
the sensitivity of the spin-chain ground state to pertur- 
bations from the system-chain coupling at these points. 
The Hahn echo is a oscillating function of time r around 
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FIG. 4: (color online) (a) Hahn echo versus time r and A with 
7=1. The other parameters chose are g = 0.1, N = 246 
sites, (b) Discrete Fourier transformation of ve(t), with the 
same parameters as in figure (a). 

A = 0. Due to the coupling to the spin-chain, the oscil- 
lation is damping, and eventually «e(t) tends to zero in 
the t — ► oo limit. The difference between cases of 7 = 
and 7 = 1 is that [n-\ t k,n^k] = for 7 = 0, but it does 
not hold for 7 = 1. This is the reason why the Hahn 
echo takes different values at these critical points. Figure 
4-(b) is a discrete Fourier transformation of ve{t) with 
the same parameters as in figure 4- (a). It would pro- 
vides us the Hahn echo in the frequency domain. The 
ground state of the XY model is really complicated with 
many different regime of behavior^), these are reflected 
in sharp changes in the Hahn echo across the line |A| = 1 
regardless of 7 (as shown in figure 5, except 7 = 0), in- 
dicating the change in the ground state of the spin-chain 
from paramagnetic phase to the others. 

Up to now, we did not consider the temperature effect. 
Finite temperature is the regime to which all experiments 
being confined, but what is the finite temperature effect 
on the Hahn echo? In the following, we shall consider this 
problem by studying the contributions of one- and two- 
particle excitations to the Hahn echo. Taking a thermal 
state Pc(0) = ^e~P Hc ((3 = u^f) as the initial state of 
the spin-chain, the Hahn echo envelope can be written as 

v%(r) = ]J\l + ^e-^{n\u{ k X k u ttk \n)\, (17) 

k n 

where \n) and f2 n denote the eigenstate and correspond- 
ing eigenvalue of H c , respectively, z is the partition func- 
tion. We shall restrict our consideration to the contribu- 
tion from one- and two-particle excitations of the chain, 
namely, 

\n) e{< fc l5(7,A)), ??l t !fei ^ fe2 | ff (7,A))}, 

for k 1 ^k 2 , or i ^ j, (18) 



2r 

1.5 




-1.5 



-1 -0.5 0.5 1 

Y 

FIG. 5: (color online) Hahn echo at time r — 20 with N = 246 
sites, and g = 0.01. 




FIG. 6: (color online) Contribution of the two-particle exci- 
tation to Hahn echo at time r = 40. The other parameters 
chosen are N = 246 sites, 7 = 0.01, and g — 0.02. 

with k, k\ and ranging from —N/2 to N/2 — 1. It is 
not difficult to show that there are no contribution from 
the one particle excitation, because X k creates or annihi- 
lates two particles with k and — k jointly. The numerical 
results presented in figure 6 show the contribution of the 
two-particle excitation to the Hahn echo, we find that the 
quantum critical points can influence the Hahn echo at a 
finite temperature. For the parameters chosen in figure 
6, the contribution form the thermal excitation is larger 
than that from quantum fluctuation when P < 72 — (3 C . 
Here we have scaled out an overall energy scale denoted 
by J. J may be taken to be of order lOOOif , that is the 
order for the antiferromagnetic exchange constant of the 
Heisenberg model. It yields T c ~ 14K corresponding to 
parameters chosen in figure 6. For the transverse Ising 
model 7 = 1, /3 C is of order 10, we obtain T c ~ 100K in 
this situation with the other parameters being the same 
as in figure 6. Notice that the study here is based on 
the Hahn echo(a dynamical quantity), this would differ 
from the investigation based on thermodynamics 1 1 5j| . We 
would like to notice that the discussion on the finite tem- 
perature effect was limited to very low temperatures, be- 
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cause only one- and two-quasiparticle excitations were 
included. Nevertheless, it is interesting because it also 
sheds light on the contribution to Hahn echo from the 
first excited states, which have the same energy as the 
ground state /5 C (0) of H c at critical points. The results 
presented in figure 6 show that those contributions tend 
to zero with T — > 0. 

In conclusion, by discussing the Hahn spin echo in the 
spin-i particle coupled to critical spin-chains, the rela- 
tion between the Hahn echo and the critical points was 
established. The relation not only provides an efficient 



theoretical tool to study quantum phase transitions, but 
also proposes a method to measure the critical points in 
experiments. Up to two-particle excitations, we have also 
studied the influence of thermal fluctuation on the Hahn 
echo, it would shed light on the low temperature (with 
respect to the overall energy scale J) effects on the Hahn 
echo. 

This work was supported by NCET of M.O.E, and NSF 
of China Project No. 10305002 and 60578014. 
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